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1. INTRODUCTION 

Nowadays, among renewable energy sources, wind power has become in the lead for producing 
electricity’s. It is playing a significant role to respond the growing demand for clean and non-toxic energy 
thanks to its multiple advantages [1]-[3]. For these reasons, the number of installations around the world has 
increased considerably. Currently, because of the electricity consumption in Algeria which should be between 
75 and 80 TWh in 2020 and between 130 to 150 TWh in 2030 on the other hand [3], the Algerian government 
provides to highlight the inexhaustible renewable resources like wind and solar energy. In this context, Algeria 
have plans to install 22.000 MW by 2030 to increase the electricity production to 40% by 2030 [4]-[6]. This 
preventive state must necessarily be achievable, because it constitutes an essential solution for Algeria so that 
the energy transition from the use of fossil fuels to renewable energies is ensured. Knowing that the exhaustion 
of fossil energies is inevitable in the near future. 

Furthermore, the wind power potential is based on several parameters [7]. So, a careful study is 
required to install a wind power conversion system. The wind estimation and the wind power evaluation has 
been carried out in many countries worldwide [8]-[15]. Bamba et al. [14] analysed the performance of a 
30 MW wind power plant in Nouakchott city, Mauritania, based on different parameters as well as the 
meteorological data analysis to give a good assessment of the available wind potential. 
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Regarding literature, Weibull model is becoming very used [16]. Hoxha et al. [10] has confirmed that 
the Weibull model is better than rayleigh model in terms of fiting wind speed variation, and this, after a 
statistical analysis of the wind speed in Kitka followed by a comparative study between the two models. After 
the study of the wind speed characteristics for Medan City, Indonesia, Suwarno et al. [17] has argued that the 
Weibull model is the better model compared to gamma and exponential distribution models. 

Thus, despite that Annaba is a coastal city which has a significant wind power potential, but no work 
has been carried out in this context, above all, that currently Algeria is strongly involved in the field of 
renewable energies. For this reason that the first part of this work aims to study the wind characteristics for 
Annaba city using real data. To determine the yearly and seasonal power wind density and to estimate weibull 
scales and shapes parameters, Weibull density function (WDF) has been used. 

Wind energy can be harnessed by wind energy conversion systems (WECS) which convert wind 
kinetic energy into electrical energy. Thanks to the remarkable advantages of variable speed WECS compared 
to those at a fixed speed (better energy capture and excellent quality) they are currently the most used. The 
remarkable development in the domain of power electronics has enabled the emergence of static converters 
(SC) in WECS. Due to its advantages: bidirectional energy transfer, small harmonic deformation of input 
currents, high power factor, lowest switching frequency and also switch voltage stress, especially attention has 
been given to the use of pulse width modulation PWM rectifier in WECS [18], [19]. Currently, to have a high 
power and sine current several SC control strategies have been developed [20]-[23]. Among these control 
techniques, the direct power control (DPC) has attracted attention and curiosity of many researchers thanks to 
its advantages [20], [24]. 

Aroussi et al. [24] has developed the maximum power point tracking MPPT control to extract 
optimum power for a variable speed wind power system as well as the DPC control to adjust the active and 
reactive power. The control technique adopted has proven a good dynamic performance and a sinusoidal 
waveform of currents and voltages under variable speed. Similarly, the same problem was treated by Beniss ef 
al. [25] has studied a DPC control strategy and the MPPT control approach in the aim to diminish the powers 
errors and the selection of the suitable voltage vector. Simulation work under MATLAB/Simulink applied to 
7.5 KW doubly fed induction generator DFIG wind system was carried out, demonstrating the efficiency of 
the used control technique. Also, in this sense, Razali et al. [26] has proposed a DPC control technique based 
on a new look-up table for three-phase grid connected to 3L-NPC. The results exposed in this work show a 
various performance like: low current harmonic distortion and unity power factor as well as a better dynamic 
response. 

The present work contains of two parts, the first is based on the evaluation of wind potential per year 
and season of Annaba city in Algeria. The assessment and analysis of wind is important in the exploitation of 
wind resources. This will unquestionably allow us to have a clear forecast of wind power potential 
characteristics, and that our results can encourage investors to establish an energy development project in the 
studied site, thus taking advantage of these metrological conditions for the installation of wind turbines. 
Analysis of measurement results shows that the wind potential importance encourages decision makers to plan 
wind conversion systems installation. The recourse to energy conversion systems for the exploitation of energy 
sources, becomes more and more required. However, knowing that the wind speed profile is random, these 
wind systems require taking into account of a control system allowing to extract the maximum power whatever 
the real climatic data. In order to ensure a good performances and unit power factor with proper power 
decoupling, the second part focuses on the application of a nonlinear direct power control strategy. The 
optimum semiconductors switching state is selected through the switching table. Indeed, the basic idea of the 
switching table revolves around on dP and dQ errors variations limited by a hysteresis band, and (8) wich is 
the position. 

The present paper proposes in section 2, an evaluation of wind energy resources and power density. 
For study the wind conversion system a modelling of the whole system is carried out in section 3. After, a 
control based on DPC is explained and programmed under MATLAB software in section 4 and section 5, 
respectively. Finally, the obtained simulation results are analyzed before ending with the conclusion 


2. ANALYSIS OF WIND ENERGY 
Because of its several advantages, Weibull distribution is the commonly used model [17], [27] for 
wind speed evaluation and analysis. The mean speed is calculated as [11], [12]: 


l yn 


Um = Elia Yi o) 


Weibull parameters: scale « c » and shape « k » can be expressed as [27], [28]: 
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Moreover, the c, k and v at the desired height are given as: 
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Wind speed data considered in this work are measured and analysed for eight (08) years (2007-2014). 
Figures 1(a) and (b) illustrate the monthly and yearly wind speed average at 20 m height. We notice that, 7.8973 
m/s is the highest value of wind speed by month in July 2008. In addition, 4.6584 m/s is the lowest value of 


wind speed by month in January 2007 and February 2008. 
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Figure 1. Wind speed by (a) months and (b) years 


2012 2013 2014 


Figure 2 shows wind speed variation by season for eight years. The highest seasonal wind speeds is 
still recorded in summer. Therefore, its value during the eight years is in the range from 4.8845 m/s in Automne 
2010 to 6.7553 m/s in Spring 2008. The seasonal and yearly characteristics of wind speeds during the eight 
years at 20m height are recapitulated in Tables 1 and 2, respectively. The maximum average speed is 
5.9638 m/s in summer, with 5.4114 m/s is the minimum average value determined in autumn. 

From Figure 3 we can deduce that the scale parameter varies from 5.3379 m/s in 2010 to 6.1321 m/s 
in 2007. Furthermore, the shape parameter is maximum in 2007, 2009 and 2011 with a peak value of 20.3664 
in 2011. Its minimal value is 4.5623 recorded in 2012. In addition, the shape parameter increases considerably 
in the middle of the year between May and October. Indeed, its maximal values are in spring and summer. 
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Figure 2. Seasonal wind speeds Figure 3. Yearly Weibull parameters 


Table 1. Seasonal wind speed and Weibull parameters at 20 m 
Period Vms M/S ks cs m/s 
Winter 5,6997 10,8408 4,0086 
Spring 5,5992 29,7804 5,6710 
Summer 5,9638 32,9582 5,8746 
Autumn 5,4114 16,9746 6,2353 


Table 2. Yearly average wind speed and Weibull parameters at 20 m 
Period vm m/s k c m/s 
2007 5,899 12,7682 6,1321 
2008 5,763 06,5449 6,1299 
2009 5,462 16,3110 5,6425 
2010 5,128 09,8721 5,3379 
2011 5,846 20,3664 6,0124 
2012 5,690 04,5623 5,9563 
2013 5,733 04,6174 6,0174 
2014 5,824 04,6526 6,0563 


3. WIND ENERGY CONVERSION SYSTEM MODEL 
Conversion systems allow the conversion of the kinetic energy into electrical energy. In the case of 
low wind speeds the gearbox increases the rotation speed in order to ensure the required power. Do to the use 


of static converters, WECS presents several simple and easy control strategies compared to other systems. 
Figure 4 show the studied WECS. 
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Figure 4. Schematic diagram of WECS 
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3.1. Mechanical part 
In (9) describes the wind power [2], [29]: 


p.s.vi 


BEGAN: (9) 


The power coefficient C, can be expressed as: 
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Where: p is the air density and R represent rayon. 
The torque is provided by the expression: 


Ptur SVP 1 
Teur = Cp (A, B).? à 
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We can model the mechanical part by : 
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When the tip speed ratio 4 is optimal Aopt and also Cyop¢. So, we can deduce: reference electromagnetic torque 


and the maximum power, determined by MPPT control [28]: 
Tem = kopt: Qa 


= 3 
a kopt On 


Cpopt-P-T- R5 

With: kopt = ET 
3.2. Modelling of the DFIG 

Considering that the grid is balanced we can write [2]: 


Pas Ls 0 Lm fias 
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Par| — Lm O Ly 0||ia 
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With [2] 
APds 
Vas =%2=0 
Vas = Ws. Pas = V 


P 2 è è è 
Te = 37 Lm (igs lar — las: igr) 


Where: 

Qs = Lm. İmsis the stator flux and @, is the rotor flux; 

L,, Lm and L,: stator, magnetizing, and rotor inductances, respectively 
V,, igs: Stator voltages and currents; 

V,, ig,: rotor voltages and currents; 

R,, R,: rotor and stator resistances; 

We, W,: synchronous and rotating angular frequencies; 

ws: synchronous angular frequency; 

Ws, = We — w: slip frequency; 
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Tę: electrical torque; 
P: number of poles 


3.3. Modelling of three-phase converter and PWM control 
The converter mathematical modelling can be written in matrix form as: 


Va La d La Va 
Va| =R [ip | + LE | io | + [Vo (22) 
Voc ic ic Ve 


The voltages V,, , Vp and V, are defined as [30]: 


Vdc 

Va = £ (2Sa — Sp — Se) 
Vac 

Vp = (Sa — 254 — Se) (23) 
Vdc 

Y= g (Sa — Sp — 2Sc) 


The semiconductors opening and closing of the inverter depend on the states of the control signals (Sy ,Sp,S¿), 
respectively, as they are defined: 


F T, closed and T; open 
Sa = 


0 T open and T; closed 


_ (1T, closed and T; open 
da lo T, open and T, closed ey 
ei E T closed and T; open 
© LOT, open and T; closed 
+1,5, =—1 
Sk = { paa 25 
K LoS = +1 (25) 
Or, k=a, b, c. 
DC-Bus equation is given by: 
adc A 4 . < 
C a Sala + Spip + Sele — Set; (26) 


4. CONTROL STRATEGY DESCRIPTION 

In order to calculate the semiconductors switching state for powers control, we used a hysteresis 
regulators. The powers variations are the inputs of the regulators. Furthermore, the basic idea of the switching 
table articulated on dP and dQ errors variations limited by a hysteresis band as shown in Figure 5. The DPC 
[20]-[23] is a nonlinear and direct control because it doesn't use any modulation technique. 
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Figure 5. Inverter control block 


Thus, there are eight combinations which develop at the inverter output, six vectors of active voltages 
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(V;with i=1:6) and two vectors of zero voltages (Vj) and (V7), as given in Table 3. The eight possible voltage 
vectors in the Clark plane (a, B) are shown by the polygon of Figure 6. They are expressed by the: 


jk 
v= ace! "a (27) 


With, k=1,..... 6 and V, =V,=0 


Table 3. Switching states and output voltages B 4 
Sa Sp Se Vea Ven Dec Uca Ucp Vi 
0 0 0 0 0 0 0 0 Vo 
0 0 1 Va Vac  —2Vac Vac Vac Vs 
3 3 3 V6 VZ 
0 1 0 Va Vac Vac Vac Vac V3 


3 3 3 J6 v2 


0 1 1 -Wa Vac Vac 2 0 Vi 
3 3 3 = | Vac 


1 0 1 Vac -2 Vac Vac Vac Vac Vo 
3 3 3 J6 V2 
1 1.0 Vac Vac —2Vac Vac Vac V2 
3 3 3 V6 V2 ; aos 
I A 1 0 0 0 0 0 Vz Figure 6. Switching polygon 


The hysteresis controllers are at two levels for the P and Q powers where their outputs are represented 
by the boolean variables dP and dQ indicating, respectively, the higher or lower overruns of the following 
active (P) and reactive (Q) power errors according to the following logic: 


Prop -P >h, >dP=1 
Prep — P < h, >dP =1 
Qref — Q 2 h¿2dQ =1 
Ore -Q < hg 3dQ =0 


(28) 


Where, h, and h, are the hysteresis bands regulators. For this purpose, in the case where the errors are greater than 
zero the regulator output is equal to 1 else its equal to 0. The purpose of this type of power controller is to maintain 
the end of the error of the active or reactive power in a hysteresis band of 2h. The hysteresis corrector output indicates 
if the power error must be increased (dP=1), (dQ=1) or decreased (dP=0), (dQ=0). Furthermore, Va can be regulated 
using PI controller as shown in Figure 7. Figure 8 presents the regulation of DC voltage with a PI controller. 


Vác ref 
Regulator 


Vac ref E(S) 
— Tm 
Prof | 


Figure 7. Calculation of the reference power Figure 8. DC voltage regulation with a PI controller 


Vac 


Where, the output voltage is multiplied by the effective DC bus voltage Vac to generate P,-¢¢. We can write [29]: 


Vac _ Kp.k¡/K 
Vacref  S?+kp/K.S+kp.k¡/K 
Where: 


(29) 
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V2CVac ref 
3Vs 


K= (30) 


According to the in (31), the relation between Vac and Va, +e reveals be a second order transfer function: 


Vdc _ wr 


Vac ref — $242E wn Stun 


(31) 


The PI controller parameters are chosen according to the following relationships: 


w 32 
ki = "zz ( ) 
The reference continued current is expressed by: 


Ted ref = Lo ref + Lon (33) 


Where: 
Ien: reference current in the capacitor given by the regulator PI of the voltage 
Ic ref: measured load current; 
Ira ref: reference continued current. 
The product of the reference DC current with the DC voltage gives the reference active power: 


Pref = Vac ref Ira ref (34) 
We can write the active and reactive powers as: 


Ps = Vea İsa F Vsp- isp + Ve ise 


1 : ; : 35 
la — z Vs» a Vo)isa + (Yc 7 Voc disp oh (Vea Von) isc ( ) 


Hence, we can estimate the voltage using the: 


l= wali x lel (36) 


The reactive power reference is maintained at the desired value (generally zero to have a unit power 
factor). The P,¿ is calculated by multiplication current value generated and the optimum voltage [31]. The two 
hysteresis controllers’ outputs and the source voltage vector position are the switching table inputs. Depending 
to the vector angle of the reference source voltage on the axis (a), we decompose the position space into twelve 
sectors as shown in Figure 9. Where, it can be determined by [32]: 


(WN — 2) 4< Oy <(N- DE (37) 


Where, N=1, 2....12. 
The sector number (N) is instantly determined by the voltage vector position given by: 


9 = arctg £) (38) 


Table 4 presents the switching table which is the essential part in the DPC control [23]. It enables the selection 
of the voltage vector to allow the orientation of the active and reactive powers towards their designated values. 


Table 4. Switching table 
dP dQ O; O» O; Os, Os Os O, Os Os O10 On On 
1 1 Vo V7 Vi Vo Va V7 Vz Vo Va V, Vs Vo 
1 0 V7 V7 Vo Vo Vz V Vo Vo Vr V, Vo Vo 
0 1 Vo V; Vi V Vo V3 V3 Va V Vs Vs Vo 
0 O Vi V V Vs V Ve Ve Vs Vs Ve Vs Vi 
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Figure 9. Position of different sectors [31] 


5. RESULTS AND DISCUSSION 

The whole system was carried out using MATLAB/Simulink (see Figure 10). We have considered a 
wind profile evolving according to the Figure 11(a) with a wind speed of 12 m.s”! between 0< t <1s and equal 
to 9 m.s”! between 1<t< 2s. Thus, the couple and the generator speed evolve according to the paces represented 
by the Figures 11(b) and 11(c), respectively. However, the electrical network voltages and currents waveform 
are purely sinusoidal whose figures are shown, respectively, by Figures 12(a) and 12(b) (see Appendix). 
Moreover, for more details, the Figures 12(c), 12(d) and 12(e) (see Appendix) present zooms or parts of these 
last figures. 

We notice that, the amplitude of the network currents changes with the change in wind speed. And 
that, since we have considered a reference of the variable reactive power (Q=0 if 0< t <0.5s; Q=4>KVar 
if 0.5< t <ls; Q=2 KVar if 1< t <1.5s and Q=0 if 1.5< t <2s) we note that the currents waveforms (Figure 12(c) 
change in amplitude (see the circled zooms A and B). Finally, Figures 13(a)-(c), highlight the paces of: active 
power, reactive power and a zoom of reactive power. We can see that the reactive power follows its reference 
perfectly as defined above. However, at times t=0.5 s and t=1.5 s, the active power keeps its reference value 
notwithstanding the reactive power variation. These, witnesses and valid the studied control strategy, which 
ensures a powers decoupling. 
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Figure 10. Schematic diagram of the proposed control algorithm 
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6. CONCLUSION 


In this paper, an evaluation of wind potential of a site in Annaba, Algeria was carried out in prevision 
of wind power production chain installation. The results showed that Annaba city has a reasonably average 
wind farm with a peak of nearly 6.7553 m/s at 20 m altitude in Spring 2008 and 7.8973 m/s in July 2008. Based 
on this analysis, it should be noted that we can use the wind capacity for power production. Finally, as a coastal 
city, Annaba can be suggested as one of the desirable’s places for the installation of wind turbines. Furthermore, 
the DPC control strategy has used. The simulation results of P and Q powers give very high performances in 
term of response time and overshoot. Furthermore, a good decoupling between the active and reactive powers 
was assured by the proposed control strategy thus guaranteeing a unity power factor for 0<t<0.5s and 
1.5< t <2s. So, this command has proven to be sufficient in such cases. However, the application of robust and 
intelligent commands will certainly be the subject of our future work. 


APPENDIX 
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Figure 12. The simulation results for, (a) voltages network, (b) voltage zoom, (c) currents network, 
(d) currents zoom for t=0.5, and (e) currents zoom for t=1 
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